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Abstract 



Screening of initial parton production due to the presence of on-shell 
partons in high-energy heavy-ion collisions is discussed. It is shown that 
the divergent cross sections in the calculation of parton production can 
be regulated self-consistently without an ad hoc cut-off, and that the 
resultant parton production and transverse energy production rate are 
finite. Consequences on the energy density estimates are discussed. 
PACS numbers:12.38.Mh,25.75+r,12.38.Bx 
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Hard and semihard parton scatterings have been shown to play an important role 
in both high energy pp(p) and ultrarelativistic heavy ion collisions. In pp(p) collisions 
they are believed to be responsible for the rapid growth of the inelastic and total 
cross sections [0,0 , and in AA collisions they are expected to produce a large amount 
of transverse energy in the central rapidity region [Q-|6|]. Because the dominant QCD 
parton cross sections are singular in the soft scattering limit, most model calculations 
of minijet production via perturbative QCD (pQCD) have introduced an infrared 
cut-off, po, corresponding to the smallest permissible transverse momentum transfer 
in a 2 — ► 2 parton scattering, to which perturbative QCD can still be applied. This 
cut-off is used in most of the models to separate perturbative hard processes from 
nonperturbative soft interactions which can be modeled, e.g., by strings or flux tubes. 
Since there is no distinct boundary between soft and hard physics, both the hard and 
soft part of the interaction in this scheme are very sensitive to the cut-off p$. 

Heavy ion collisions differ from pp collisions in that minijets are produced in large 
number, so that a medium of minijets is formed. In the space-time evolution of a 
heavy-ion collision soft particle production will be completed after semihard and hard 
processes. Therefore, the soft interactions are expected to be screened by interactions 
with the semihard quanta (minijets). We here take a step further by proposing that 
the perturbative semihard particle production should also be screened by the processes 
which have happened even earlier, i.e., by the harder processes. In this paper we will 
study the influence of the final-state interactions among minijets. We will investigate 
to what extent minijet production is screened and whether the screening effect can 
provide an effective self-consistent regularization of the infrared behavior of minijet 
production. 

Consider an example of two hard processes as illustrated in Fig. 1. Let us assume 
jets from the first hard scattering are produced in the 90 degree direction and have 
a large transverse momentum p-ri- Because of this large transverse momentum, the 
interaction point is well localized in the transverse direction to a distance of 1/pri- 
As these two jets travel in the transverse direction, they will experience secondary 
interactions. These secondary interactions can give rise to many nuclear effects of hard 
scatterings, i.e., energy loss and Cronin effects. Since we are interested in the screening 
effects caused by these final state interactions, we will consider the interactions of the 
produced hard partons with the propagating parton in another scattering nearby as 
shown in Fig. 1. A semiclassical estimate of the screening requires that different 
scattering events can be treated as incoherent. We will show that this condition is 
satisfied if the produced partons, which screen other softer interactions, can be treated 
as on-shell particles. 

In a parton model, we treat partons from the nuclei and out-going partons as on- 
shell. Therefore, the propagator in the first hard scattering must be off-shell by —p\\- 
One notices that there is one extra loop in the Feyman diagram in Fig. 1 as compared 
to that of two independent scatterings. The loop-momentum integration will then 
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put one of the propagators on shell. Since we demand that the second scattering also 
has transverse momentum transfer p?2, the propagators in this scattering will also 
be off-shell by — p\ 2 - Therefore, the dominant contribution to the loop-momentum 
integration comes from the pole of the propagator between the two scatterings in 
which the propagating parton will be put on-shell. The condition for this contribution 
to be dominant is that the distance Axj_ between the two scatterings must be larger 
than the interaction range of the two hard scatterings which are determined by the off- 
shellness of the exchanged gluons, i.e., Ax± > max(l/pn, Vpt2)- ^ this condition is 
not satisfied, the propagating parton between two scatterings cannot be treated as real 
and consequently one cannot treat the multiple scatterings as incoherent scatterings. 
One can thus consider 



as the formation time of the produced parton in the mid-rapidity from the hard or 
semihard scatterings after which they can be treated as real (on-shell) particles and 
can screen other interactions with smaller transverse momentum transfer. Here we 
used a dimensionless coefficient a of the order unity to characterize the uncertainty 
in the formation time. 

Let us visualize high-energy AA collisions as in Fig. ^. In this idealized space-time 
picture the incoming nuclei are very thin slabs, the primary semihard parton-parton 
collisions all start at t = and the forming system at central rapidity is longitudinally 
boost invariant. In the space-time evolution of the collision the partons with larger 
Pt are produced earlier, as implied by Eq. [l]. These hard partons will then screen 
production of partons with smaller p? later in time and space. For the central region 
around y = 0, we will consider the screening effect of partons within a unit rapidity 
interval, Ay = Y = 1. We assume minijets have a plateau in the rapidity distribution, 
which is a good approximation both at LHC and RHIC energies for our choice of Y. 
For a fixed px, partons with larger rapidities form later in our given frame, along the 
proper time hyperbola as in Fig. [|. In this picture, the formed partons can therefore 
only screen those scatterings which are in the same rapidity region. 

Following Ref. JFj, we can estimate the static electric screening mass generated 
by the produced minijets. The number distribution of minijets produced in an AA 
collision at an impact parameter b = can be written as |4| 






(2) 



where Taa{^>) is the nuclear overlap function and 




(3) 
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is the minijet cross section. The hats refer to the partonic sub-processes, %i is the 
momentum fraction of the initial state parton i, is the transverse momentum, 
yk is the rapidity of the final state parton k, fi is the parton distribution function, 
and K = 2 is a factor simulating the contribution from the 0(a^) terms in the 
cross section In what follows, we will compute the minijet production from the 

complete lowest order formula in Eq. |3| but will finally treat all the minijets as gluons. 
This should again be a good approximation, since gluons clearly dominate the minijet 
production at these energies and transverse momenta 

To obtain an estimate of the average parton number density in the central region at 
a given time Tf, we divide Eq. ||by the approximate volume V = nR\Az ~ 7r -R^t/ Ay 
of the produced system. Then the screening mass from Ref. Jj]] becomes 

" B^tv 2 -*h(V72)] 1^ 1^l o (4) 

At this point there are several remarks to be made. First, as done in Ref. we 
have replaced the thermal (isotropic) Bose-Einstein distribution by the non-isotropic 
minijet distribution, and assumed that the screening mass can still be calculated from 
the equilibrium formula given in [0]. 

Second, note that the lower limit in the /^-integration in Eq. f| is pt- In this way, 
we are considering a maximal screening effect by assuming that all the quanta with 
transverse momenta > Pt = a / T screen the formation of partons at transverse 
momenta kx < Pt- Also, as explained in the previous paragraphs, only the quanta 
within the rapidity window Y are assumed to contribute to the screening mass. How- 
ever, from Eq. |4] we see that as far as the interval Y is small enough, the screening 
mass is almost independent of the choice for Y, due to the division by the volume of 
the system. Finally, the scale choice for the running coupling is not clear. Here we 
just use the pt of the primary semihard processes. 

As we know from studies of pQCD at finite temperature, perturbation theory 
does not generate static screening in the magnetic part of interactions |TI|]I2]| . One 
normally has to introduce a nonperturbative magnetic mass to regulate the infrared 
singularities in the transverse part of a gluon propagator. Since the magnetic screen- 
ing mass in the static limit is not known, let us proceed in the following phenomeno- 
logical way to get an order of magnitude estimate on the screening effects. We will 
use the computed electric screening mass as a regulator for the divergent i- and u- 
channel sub-processes. We will simply make a replacement —i{u) — > —t(u) + n 2 D in 
the minijet cross section of Eq. |[ Since the s-channel processes are not important 
here, we do not need to make the above replacement for s. Then we are also safe 
from the contradiction that the invariant s is always positive but the screening mass 
Hd is computed for a spacelike gluon. 

In Fig. [3| we show the screening mass \id and the screened one-jet cross section as 
functions of pt- In the upper panel the results are shown for the LHC energy (per 
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nucleoli pair) y^s = 5.5 TeV and in the lower panel for the RHIC energy yfs = 200 
GeV. We start our numerical computation at sufficiently high pt (i.e. at early times 
Tf) where the screening mass is small enough and the cross section is not screened in 
any significant way. At this Pt we then compute the screening mass [a%(pt) from Eq. 
d and with the obtained (i d (pt) we compute the screened cross section at pr — Apr 
(with a sufficiently small step Apy). From the resultant cross section we get the 
new screening mass at px — Apr and so on. This iterative procedure gives us a self- 
consistent picture of screening: the bigger the cross section becomes the denser the 
parton medium gets and the larger a screening mass is generated, slowing down the 
rise of the cross section towards smaller px- In this way, the medium of produced 
minijets regulates the rapid growth of the jet cross section as seen in Fig. Finally, 
at {Id ~ Pt, the cross section saturates. From Fig. ||] we see that for sufficiently 
large nuclei and for high collision energies this happens in the perturbative region 
Pt Aqcd- To study the sensitivity of the results to the uncertainty relation (1), we 
show curves corresponding to a = 1 and a = 2. Note that due to the self-consistent 
iteration procedure the resulting difference in ht> is not a trivial multiplicative factor 
V2. 

In the computation we have used the MRSA parton distributions fi(x, Q) with a 
scale choice Q = Pt H3 . At this point it is worthwhile to note that the singular parton 
distributions are stable against scale evolution, and the small-x growth xf\ ~ x~ s 
persists even when going down to scales pt < 2 GeV. Most of the cross section of 
Eq. |3] at fixed pt comes from momentum fractions x ~ 2pt/\^s, so that a smaller 
scale pt also corresponds to a smaller x. Therefore, when going towards smaller pt, 
there is a competition between the scale evolution (causing a decrease in the parton 
density) and the small- x rise. As a result, the parton luminosity does not decrease 
fast enough to cancel the divergence in the t- and u-channels at small pt- This is 
basically the reason for the rapid growth of the non-screened cross sections as seen in 
Fig. i 

In the top panel of Fig. [| we show the average number of produced minijets with 
Pt > Po an d \y\ < 0.5 in a central AA collision with A = 200. The curves are shown 
as functions of po and they are obtained simply by integrating the cross sections in 
Fig. |3| over pt from p upwards and multiplying by T^(0). The average transverse 
energy carried by these partons is obtained by weighting the cross sections of Fig. |3] 
by pt and integrating again over pt in the region p < pt < y/s/2. The corresponding 
curves are shown in the lower panel of Fig. |] as functions of Pq. Note also that po 
here gives the time for formation of the system according to Eq. [I]. 

We want to point out how the parameter p is needed as a cut-off only for the 
unscreened cross section. This cut-off makes the unscreened minijet production finite 
but the results depend strongly on the cut-off, as seen in Fig. [4[ When the self- 
consistent screening is included the cross section saturates and there is no need for 
such a cut-off anymore, provided that the saturation occurs in the perturbative region 
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as it does in Fig. |[ However, one may still ask how many minijets are produced into 
the central rapidity unit with transverse momentum larger than some lower limit po, 
and what is the average transverse energy from these quanta. Fig. [| gives an answer 
to these questions. Especially at the LHC energies, where the number of produced 
minijets becomes very large, and the parton system very dense, the average transverse 
energy practically saturates near p ~ 1 GeV. 

The weaker po-dependence should also make the perturbative energy density es- 
timates of the early minijet system at T/(po) = a /Po more reliable. In Fig. ^| the 
transverse energy density |TJJ 
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is plotted against p$. Also the corresponding formation time is shown simultaneously. 
It is interesting to notice how the maximum energy density is reached at earlier times, 
and how the expansion of the system starts to dominate the evolution earlier when 
screening is taken into account. One should, however, keep in mind that the finite 
overlap of the colliding nuclei will reduce the energy density estimates especially at 
RHIC energies |l5j . On the other hand, roughly half of the released final transverse 
energy is expected to originate from the soft particle production at RHIC energies 
PJ5|JT5[| . This will slow down the decrease in e at a later stage [15[ . 



The screening effect we discussed here is due to final state interactions of produced 
partons. We should point out that there will also be screening effects, coming from 
initial state interactions or the color confining radius in the dense parton clouds 
of a high energy nucleus fTE] . Ideally, one should take into account both screening 



effects in the calculation of parton production. Nuclear shadowing due to initial-state 
interactions will also reduce the estimates on semihard parton production especially 
at the LHC nuclear collisions, perhaps even at RHIC energies P,|T7|. Nuclear gluon 
distributions are still uncertain to a large extent, and especially the scale dependence 
of nuclear gluon shadowing has to be studied carefully, e.g., as done in JC7J, but with 



the modern parton distributions like used in this paper [|T^ . 

In summary, we have considered color screening of initial semihard parton pro- 
duction in a phenomenological but self-consistent manner. Following Ref. [0] we have 
computed the static electric screening mass of the parton system by using the number 
distributions of produced minijets and by taking into account the formation time of 
the minijets. We do not try to estimate the corresponding magnetic screening mass, 
instead we use the obtained electric screening mass as a regulator for the divergent t 
and w-channel sub-processes. We have iteratively computed the lowest order minijet 
cross section with a feedback from the screening mass towards the region p^ ~ Hd ~ 2 
GeV. We have shown that for sufficiently large nuclei and for sufficiently large en- 
ergies, we obtain saturation of the minijet cross section in the perturbative region 
Pt ^> h-QCD- From this we conclude that the average number of minijets in the 
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central rapidity unit is finite and, furthermore, there is no need for an ad hoc cut-off 
parameter p . However, the question of how many minijets are produced with p? > Po 
is still quite valid, and we have shown that the answer depends on the lower limit 
Po much more weakly than when screening effects are neglected. In the estimates of 
perturbative transverse-energy production there is practically no po-dependence be- 
low p < I GeV at the LHC energy. Finally, we have shown how much the screening 
reduces the Bjorken estimate of the energy density of the early parton plasma. 
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FIGURES 

FIG. 1. An example of a process contributing to the screening. Off-shell gluons arc 
drawn in thick lines. 

FIG. 2. Longitudinally boost invariant space-time picture of an AA collision. The 
z-axis is along the beam, and the longitudinal proper time is r = \/t 2 — z 2 . Partons in the 
central region within a rapidity interval Y (between the solid lines) and with transverse 
momenta p? > a/r\ have formed earlier than n, and thus screen the formation of partons 
at r > r\. 

FIG. 3. (a) Differential minijet cross section da^/dprdy at y = and screening mass 
Hd as functions of transverse momentum px in a pp sub-process of a central A A collision 
at y/s = S.S^TeV with A = 200. The dashed line is for unscreened cross section, the solid 
and dotted curves for the screened case with different parameters a from Eq. 1. Nuclear 
shadowing is not taken into account, (b) The same as in panel (a) but for y/s = 200^4GeV. 

FIG. 4. (a) The average number of minijets produced with px > Po and \y\ < 0.5 in a 
central AA collision as a function of po (see the text). The solid, dotted and dashed curves 
correspond to the curves in Fig. 3, and the two sets of curves are for energies y/s = 200 and 
5500 ^4GeV. (b) The average transverse energy of the minijets of panel (a) as a function of 
Po (see the text). Labeling of the curves is the same as in panel (a). 

FIG. 5. Transverse energy density e of the perturbatively produced minijet system 
as a function of the smallest transverse momentum p$. The corresponding formation time 
axis is also shown. The upper set of curves is for yfs = 5500 AGeV and the lower set for 
y/s = 200 ^4GeV. The solid and dashed lines show the estimate with and without screening, 
respectively, for a = 1. 
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